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Abstract 



Using our solution for the open spin- 1/2 XXZ quantum spin chain with N spins 
and two arbitrary boundary parameters at roots of unity, the central charge and the 
conformal dimensions for bulk hole excitations are derived from the correction to 
the energy (Casimir energy). 



1 Introduction 



The integrable open spin-1/2 XXZ chain has been subjected to intensive studies due to its 
growing apphcations in various fields of physics, e.g., statistical mechanics, string theory 
and condensed matter physics. However, obtaining exact solutions for this model has been a 
rather challenging and elusive task for many years. Various progress have been made in ob- 
taining solutions for this model, either using the Bethe ansatz approach for diagonal PP-[1], 
constrained nondiagonal [5]-[9] and nondiagonal cases at roots of unity [I0]-[T3], or using 
the representation theory of the g-Onsager algebra for general nondiagonal cases [H]. Ap- 
proaches based on boundary Temperley-Lieb algebra and its representations have also been 
presented recently, from which the spectral properties of the chain have been studied |15j . 
Upon obtaining the desired solution, the next natural question that needed to be addresed 
is its practicality within various contexts. One important area where these solutions have 
found creditable applications is in determining finite size corrections to the ground state 
energy. By relating to conformal invariance, these finite size corrections are shown to be 
related directly to other crucial parameters like the critical indices, central charge and con- 
formal dimensions [I7j-[20]. There are few methods and approaches to accomplish this task. 
De Vega and Woynarowich [21] derived integral equations for calculating leading finite-size 
corrections for models solvable by Bethe Ansatz approach [22]. This was then generalized to 
nested Bethe Ansatz models as well [23j. Another approach was introduced by Woynarowich 
and Eckel [211 [25] , which utilizes Euler-Maclaurin formula and Wiener-Hopf integration to 
compute these corrections for the closed XXZ chain. Others have also studied more general 
integrable spin chain models e.g., XXZ diagonal [21 13], nondiagonal cases [2H], quantum spin 
1/2 chains with non-nearest-neighbour short-range interaction [26] and XXZ(l/2, 1) which 
contains alternating spins of 1/2 and 1 [27J, within similar framework. Other approaches e.g., 
based on NLIE (Nonlinear Integral Equations) have also been successful in determining these 
effects for integrable lattice models [29] and related integrable quantum field theories, such 
as the sine-Gordon model with periodic [30]-[32], Dirichlet [33]-[37] and Neumann boundary 
conditions [28lli2]. 

With similar aim in mind, utilizing an exact solution for the integrable spin-1/2 XXZ 
chain with nondiagonal boundary terms at roots of unity we found earlier for even number 
of sites [m [38] 0, and extending the solution to account for odd number of sites as well, we 
compute the correction of order 1/A^ (Casimir energy) to the ground state energy together 
with its low lying excited states (multi-hole states). We employ the method introduced 
by Woynarovich and Eckle [24j that makes use of Euler-Maclaurin formula [15] and Wiener- 
Hopf integration [36]. In particular, we compute the analytical expressions for central charge 

^This solution, in contrast to [B|-[5] does not assume any constraint among the boundary parameters. 
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and the conformal dimensions of low lying excited states. We also compare these analytical 
results to corresponding numerical results obtained by solving the model numerically for 
some large number of sites. 

The outline of this article is as follows. In Section 2, we review the Bethe Ansatz solu- 
tion [HI [38]. We also present an extension of that result to include solution for odd A^. In 
Section 3, we present the calculation of correction to the ground state energy and hence 
our results for the central charge and conformal dimensions of low lying excited states. We 
notice that the lowest energy state for even N of this model has one hole. Hence, the true 
ground state (lowest energy state without holes) lies in the odd sector. Similar behaviour 
are also found for the open chain with diagonal boundary terms, for certain values of bound- 
ary parameters |47j. It is known that (critical) XXZ model with nondiagonal boundary 
terms corresponds to (conformally invariant) free Boson with Neumann boundary condition 
whereas the diagonal ones are related to the Dirichlet case [311 ESI ESI |12]. Although the 
model we study here has nondiagonal boundary terms, we find that the conformal dimen- 
sions for this model resemble that of the Dirichlet boundary condition. Some numerical 
results are presented in Section 4 to confirm and support the analytical results derived in 
Section 3. Here, we solve the model numerically for some large but finite A^ and further 
employ an algorithm due to Vanden Broeck and Schwartz [39]-[l0] to extrapolate the results 
for A^ ^ cxD limit. We conclude with a discussion of our results and some potential open 
problems in Section 5. 

2 Bethe Ansatz 

We begin this section by reviewing recently proposed Bethe Ansatz solution [TTl [38] for the 
following model [13|, [B] 

where the "bulk" Hamiltonian is given by 

^0 = 9 5Z + ^n^n+l + cosh f] . (2.2) 

^ n=l 

In the above expressions, a^, a^, are the usual Pauli matrices, rj is the bulk anisotropy 
parameter (taking values t] = with p odd), a± are the boundary parameters, and A^ is 
the number of spins/sites. Note that, this model has only two boundary parameters. Other 
boundary parameters (as they appear in the original Hamiltonian in [13]) have been set to 
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zero. We restrict the values of a± to be pure imaginary to ensure the Hermiticity of the 
Hamiltonian. The Bethe Ansatz equations for both odd and even are given by 



J = 1 , 2 , . . . , Ml , 
j = l,2,... ,M2. (2.3) 



where 



r/ \ 4/ , , / \\2N sinh2u sinh(2M + 4r7) 

oiu) = 2 (smhusmhfu + 2r7)) — — r — smhf-u + 77 + a_) 

sinh(2M + 1]) sinh(2M + Sr]) 

sinh(M + 7] — q;_) sinh(M + r] + a^) sinh(M + r] — a^) cosh'^(M + 77) (2.4) 



and 



Ma 

Qa{u) = Ylsmh{u-uf'')smh{u + uf^ +ri) , a = 1,2, (2.5) 
i=i 

Ml and M2 are the number of Bethe roots, u^^^ and u^^^ (zeros of Qi{u) and Q2{u) respec- 
tively). However, h^^\u) and h^'^\u) differ for odd and even values of A^, as will be noted 
below. The energy eigenvalues in terms of the "shifted" Bethe roots m^"'' are given by 

E = -sinh%VV ^- + -(A^ - l)coshr7. (2.6) 

2 'tt^^sinh(nf)-2)sinh(nf) + 2) 2^ ^ ' 

where u['^^ = + ^ 



2.1 Even N 

We begin by recalling |3B] the structure of roots distribution for this case. The Bethe roots 
u^^^ for the lowest energy state have the form 

/iAf'^^ : J = l,2,...,M(,i) 

where Xf'^^ are real. Here, M(i,i) = M(2,i) = f , and M(i,2) = 2±i, M(2,2) = The fiXf'^^ 
are the zeros of Qa{u) that form real sea ("sea roots") and /iA[,"'^^ are real parts of the "extra 
roots" (also zeros of Qa{u)) which are not part of the "seas". Hence, there are two "seas" of 
real roots. We employ notations similar to the one used in [28] . 

M,(,(A+¥)) c°^h(MA + f)) 

sinh (/i(A - f )) ^ 2/x^ cosh (/i(A - f )) ^ ' 
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Rewriting bulk and boundary parameters [28], 7] = i/j,, a± = i/iaj§, where /i = and 
taking 



p+i 



8sinh2^+i {u + 2ri) cosh^(u + rj) cosh(u + 2ri) 



sinh(2-u + 3r]) 

the Bethe Ansatz equations (12.31) for the sea roots then take the following form [11], |38] 
e,(A»''')^"+'L(A»''')c,+2._(A;"')e,_2._(Aj"')e,«„,(Af'')e,_,„,(A;">)l"' (2.10) 



n h(M"' - ").2(a;»' + Af 



fc=l 



(p-l)/2 



n [*(Af"-Af=')*(A5"'+At'^'') 



fc=l 



and 
ei(A 



Af/2 



(2.11) 



(P+l)/2 

.1 n *(Af"-Ai''=')*(Af"+A^^') 



k=l 



= -nh(Ar^-A^^))e,(Af'nA(^'^) 
fc=i 

respectively, where j = 1 , . . . , The corresponding ground-state counting functions are 
h^'^A) = ^{(2iV + l)gi(A) -n(A) -gi+2a_(A) -gi_2a_(A) -gi+2a+(A) -gi-2a+(A) 



N/2 



k=l 



5]|*(A-Af») + *(A + A 



(P~l)/2 



Ar) + r.(A + AP; 



fc=l 



}(?-12) 



and 



h(2)(A) = i^^(2iV + l)gi(A)-ri(A) 

Af/2 



A;=l 



(P+l)/2 



5^ g2(A-A(^'^)) + g2(A + Ai^'^)) - r,(A - Ai^'^)) + r,(A + A 



fc=i 



2.13) 



where gn(A) and r„(A) are odd functions defined by 

9n(A) = TT + z lne„(A) = 2 tan~^ (cot(n/i/2) tanh(/iA)) 

r„(A) = i\ngn{\). 
These counting functions satisfy the following 



(2.14) 



(2.15) 



In f l2J[5ll above, / = 1 , 2 



^ Bethe Ansatz equations written in this and subsequent sections are true only for suitable values of a±, 
< \a±\ < ^ 



namely < \a±] < , a+a_ > , where i/ = p + 1 



2.2 Odd N 



In this section, we present an extension of the previous results to include solutions for 
odd values. The roots distribution is similar to the previous case, but now we have 
= M(2,i) = and M(i,2) = M(2,2) = Using the following in Q, 



u 



sinh(-u — a+ + 77) sinh('U + a+ + 77) sinh^^"'""'^(M + 2?]) cosh^(M + rj) cosh('U + 21]) 



sinh(2M + 3?7) 



we obtain the Bethe Ansatz equations 



(2.16) 



.(1,1) 



(l,l)^ 



.(1,1)- 



(2.17) 



(N+l)/2 (p-l)/2 

n h(Ay'"-Af")e,(A<"'+Af")l n [^.(Aj"' - Ar'te(Af ^" + A 



k=l 



k=l 



^(1,1) ,,(2,2)^_^^(1,1) , ^^2,2)^ 



and 



ei ( Af V^'-'M ^71 ( Ar "Oei+2a, ( A}-0ei-2a, ( A;-0 



.(2,1) 



(2,1) ^ 



,(2,1)- 



(2.18) 



(Af+l)/2 (p-l)/2 

n h(A;"'-Ai"')e,(Af" + Ai"')l n U(Af " - Ai'^^'te(Af " + A 



k=l 



k=l 



^(2,1) ^(1,2)^_^^(2,1) , ^a,2)^ 



respectively, where j = 1 , . . . , Note the presence of parameter-dependant terms in 

both the equations above. One can also notice the number of extra roots changes from 
to 2^ for Qi{u). The ground-state counting functions for this case read 

h^'^A) = i-|(2Ar+ l)gi(A) -ri(A) -gi+2a_(A) -gi_2a_(A) 



(Af+l)/2 



and 



k=l 



1 



(p-l)/2 



- Yl [?2(A-Af^)) + g2(A + Af^))]- Yl [r2(A-Al^'^))+r2(A + Af^))](}.,19) 



k=i 



h^'HA) = —[i2N+ l)gi(A) - n{X) - gi+2aJA) - gi-2a+(A) 



(Af+l)/2 



- E ^A-Al^'^))+g2(A + Al^'^) 

k=l 

As for even A^, we again have the following 

h«(A,)=J, J = l 
5 



(p-l)/2 

Y h(A-A(^'^))+r2(A + A 



k=l 



(l,2)^ 
k I 



•?0) 



Ar+ 1 



(2.21) 



where 1 = 1,2. Note that fl2.15p and fl2.2ip can be written more compactly as 

h^'HA,)=j, J = 1,...,L^J (2.22) 

where [...J denotes the integer part and fi\^N+i^ is the largest sea root for that "sea". 
Subsequently, we shall denote largest sea roots as /lAj. 



3 Finite-size correction of order 1/N 

In this section, we shall compute the finite-size correction for the ground state and low lying 
excited states. For these excited states, we restrict our analysis to excitations by holes which 
are located to the right of the real sea roots. Applying (12.71) to (12. 6p . we get the lowest state 
energy eigenvalues for chain of finite length N, 

N+l 



+ ^{N -l)cosfi. (3.1) 



where notations from [2H] have again been adopted 

Id /i sin(n/i) 

«n(A) = 7^7Tgn(A) - 



27T dX TT cosh(2/iA) — cos(n/i) ' 

, Id u sin(nu) 

^ = ^;7Y^" ^ =--— fT^rrfr^^^- 3.2 

2tt d\ IT cosh(2/iA) + cos{n^) 

Note that M(a_2) in (13. ip . refers to number of extra roots for Qa{u). The first and third 
terms in the curly bracket of (13. ip are summed over the number of sea roots and extra 
roots respectively. As one considers next lowest excited state, the number of sea roots and 
extra roots change. Hence, for these states of low lying excitations (with real sea), the very 
same term in the first sum will again be summed over accordingly between approriate limits 
dictated by the number of sea roots. As for the summation over extra roots, the function 
summed over depends on the imaginary part of these roots, especially in the presence of 
2-strings. However, as one shall see, for correction (in the N ^ oc limit), only the sum 
over the sea roots contributes. The second sum in (13.111 contributes to order 1 correction 
(boundary energy) which we have considered elsewhere 



[38]. 



^Equation [4.26] for the boundary energy in [55] holds both for even and odd values of N 
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3.1 Sum-rule and hole-excitations 



Now we present some results based on the solution of the model fl2.ip for = 2 , 3 , . . . ,7. 
We begin with even N case. We find for even N, excited states contain odd number of holes 
for each Qa{u). This can be seen from the following analysis on counting functions. For 
the lowest energy state the counting functions are given by fl2.12p and (12.131) . By using the 
fact that qn{X) — > sgn(n)7r — jin and r„(A) —fin as A ^ cxd and p*-'-* = -^^^ we have the 
following sum rule 

d\p">(\) = -(h<'>(oo)-hm(A,)) 

fiAi refers to the largest sea root. As before I = 1,2. We make use of the fact that 

urn/ N N 3 
^^22 

h^'^AO = f (3.4) 

From (13.31) and (13. 4p . we see that there is one hole located to the right of the largest sea 
root. Similar analysis for low lying (multi-hole) excited states yields the following 

rfA/)(A) = _(h«(oo)-h«(A,)) 

= Jjil + N^) (3-5) 

where is the number of holes (odd) to the right of the corresponding largest sea root. 
To illustrate the results above, we consider the following low lying excited states with y — 1 
and ^ — 2 sea roots and therefore different number of extra roots than the lowest energy 
state! The former case is found to have o,re hole with ^ and ^ extra roots i,r add.tion 
to a 2-string from each of the Qi{u) and Q2{u) respectively. From, 

^^22 

h«(A,) = |-l (3.6) 

one has 

l(h«(oo)-h(')(AO) = ^(^ + l) (3.7) 



''The lowest energy state has ^ sea roots. As for the extra roots, there are and of them for 
Qi{u) and (52 (u) respectively 
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Hence giving = 1. The later case has three holes with and extra roots and a 
2-string from each of the Qa{u) with a = 1 ,2. Similar analysis, 

h«(A,) = |-2 (3.8) 



yields 



l(h«(oo)-h(^)(A0) = ^(^ + 3) (3.9) 



giving Nh = 3. The total number of roots are the same for all these states. There are also 
excited states with equal number of sea and extra roots as for the state of lowest energy, 
but with position of the single hole nearer to the origin than that of the lowest energy state, 
suggesting the usual bulk hole-excitation scenario, Ehoiei^^"'^) increases as A*^"^^ — > where 
Ehoiei^^""^) is the energy due to the presence of holes and \^"'\ with a = 1,2 denote the 
positions of the holes in both "seas". We shall compute the explicit expression for energy 
due to holes shortly. 

As for the odd A^ case, we have the true ground state, namely state of lowest energy 
without hole. From the counting functions, (12.191) and (I2.20p . we have 



(3.10) 



' 1 

rfAp«(A) = _(h«(oo)-h«(AO) 

1 

2iV 

As before / = 1 , 2, and we make use of the fact that 

h«(oo) = f + 1 

h«(AO = ^ (3.11) 

From (13. lip , we see that this state of lowest energy for odd A^ has no hole, signifying the 
true ground state. Similar analysis for low lying excited states yields the following 

rfAp«(A) = ^(h«(oo)-h«(AO) 

= ^(^ + ^^) (3-12) 

where A'^^^ is the number of holes (even) to the right of sea roots. Hence, for odd A^ case, 
there are even number of holes (for each Qa{u)), with a = 1 , 2, for the excited states, e.g.. 



8 



for the first excited state with sea roots 



h^'Voo) = h - 

^ ' 2 2 

hW(A,) = ^ (3.13) 

which signifies the presence of two holes. 

It is known for simpler models of spin chains e.g., closed XXZ chain that even number 
of holes are present in chains with even number of spins and vice versa. Hence, the true 
ground state (lowest energy state with no holes) for these models is found to lie in even N 
sector. The reverse scenario (one hole in the lowest energy state for even and ground 
state in odd N sector) we find here for this model can be explained using some heuristic 
arguments based on spin and magnetic fields at the two boundaries, similar to the one given 
in Section 3 of Q In footnote 2, we notice the signs of a+ and a_ must be the same 
for boundary parameter region of interest. Hence, in Hamiltonian fl2.ip . the direction of the 
magnetic fields at the two boundaries are also the same (Both up or both down). This upsets 
the antiferromagnetic spin arrangement at the boundaries, favouring spin allignments along 
the same direction at the boundaries for chains with even N. This causes the following: 
presence of odd behaviours in the even N chain, namely the lowest energy state for even 
N sector has one hole for each Qa{u). Spins at the boundaries for the odd N chain will not 
experience such spin upset since the parallel magnetic fields favours the antiferromagnetic 
arrangement of an odd chain. Therefore, the lowest energy state for odd A^ chain has no 
holes. In other words, the true ground state exists in odd A^ sector. Further effects are the 
presence of odd and even number of holes in chains with even and odd A^ respectively as 
shown in the analysis above. 

Now, the energy due to hole excitations can be presented. We consider first the lowest 
energy state for even A^ case with one hole. Using 

N_ 

^ E ^(^ - ^i"''^) ^ / /nA')^7(A - A') - -g{\ - A^'^)) (3.14) 
for some arbitrary function g[\) and 

where I = 1,2, yuA^"'^"* = sea roots, with a = 1,2, and /iA^"^ = position of the hole for each 
of the Qa{u), one can write down the sum of the two densities 

POO 

p«(A)+p(2)(A) = 4ai(A)-/ rfA' (p«(A0+p('nA'))«2(A-A0 



^Readers are urged to refer to Figures 2 and 3 in that Section 
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+ ^[a2(A - A«) + a^iX - A^^))] + l[2ai(A) + 2a2(A) - 2b,{X) 

— Ctl+2a_ (A) — ai-2a_ (A) — ai+2a+ (A) — ai-2a+ (A) 

£-1 

- X](62(A-Al^'^)) + 62(A + Ap))) 



A;=l 

p+1 



- ^(62(A-Ai^'^)) + 62(A + Ai^'^)))] (3.16) 
Defining ptotai (A) = p(^)(A) + p^^\\) and solving (l3l6D using Fourier transform n, we liave 

Ptotali^j) = 4:S{uj) + ^R{uj) 

+ ^JH(e*"^*'' +e'"^'") (3.17) 

wliere ptotaii^^) , cl2{uj) and s{uj) are tlie Fourier transforms of ptotai(\) ,a2(A) and i^^^^^x) 
respectively. Also J{uj) = i • -R(i^) is the contribution from the second square bracket 
in (I3.16p . which will not enter the calculation for EhoieiX^""^) and will be omitted henceforth. 
The Fourier transform of hole density are the third and the fourth terms in (13.171) . which 



gives 



PhoUX) = ^[J(A - A«) + J(A - A(2))] (3.18) 

Using approximation ( 13.14^ in ( 13. ip . and making use of ( ]3.18p . one has 

^ A^TTsinu f°° 

EhoUX^"^) = / d\ a,{\)phoU\) 

■^P- J -oo 



2 



a=l 

which after some manipulation yields 



+ "-^j:-^m (3.19) 



2 

EhoUX^^) = —r-^} ^ 3.20 

^ ' 4/i ^coshTrAW ^ ' 
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Our conventions are 



/(c.) ^ / e-^ /(A) dA , /(A) - ^ / e-'"^ f{uj) dco . 



, . , , sinh ((i/ — |n|)ci;/2) , , 

a„H = sgn n , ' ^/ \ < n < 2^.. 

smn [uuj/2) 
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Generalizing the derivation to a number of holes, one has 

P.oaA) = ^EE'^(^-^i"^) (3-21) 

a a=l 

and finally the following for the energy 



= (3.22) 

„ coshTrA^^ 

Note that EhoieC^^a^) increases as A^"'' — > as mentioned above in paragraph following (13. 9p . 
3.2 Casimir energy 

In this section, we give the derivation of correction (Casimir energy) to the lowest energy 
state, for the even N case (with one hole). This result is then generalized to include odd 
N values as well as the low lying (multi-hole) excited states. We begin by presenting the 
expression for the density difference between chain of finite length (with N spins), p^^(A) + 
pj^^(A) and that of infinite length, poo(A) 

JV 

/oo 1 2 

-OO ^ AT 



JV 



/oo -| 2 

'OO ^ AT- 



as (A - 7)[p!;^(7) + P%\l) - Poo (7)] (3.23) 

In f l3.23p and henceforth, only terms that are crucial to the computation of correction 
are given. Other parameter dependant terms that contribute to order 1 correction have been 
omitted here §. Solving (13.231) yields 

JV 

/oo 1 2 

dip{x-i)[- J2 %-A?''^)-p!i^(7)] 

JV 

/oo -| 2 

rf7P(A-7)[T7 E %-Ar^)-p!?^(7)] (3.24) 
-oo „ ^f 



13- ^ 



'See |38] for details 
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where poo(A) = = 4s(A) and p{X) = ^ j^^du e ''^^YT^fe Similar equation ex- 

pressing the energy difference between finite and infinite system is also needed to compute 
Casimir energy. This is given by 

JV 

E.~E^ = -^^{ r<iAa.(A)|l f: i(A-A!,"')-pS'(A)] 

P~ 2 

JV 

/OO -| 2 

/CO 

dX a^{\)[p%\\) + p%\\) - p^{\)]\ (3.25) 
-OO 

Using (13.241) and the fact that p{ijj)ai{uj) = s{uj)a2{uj), we have 

En-E^ = -^^^^{/_^rfA5«(A)p«(A)+ (3.26) 

where S^i\X) = ^ EJ^.n 5(A - 4''^) - p!i^(A) and pS(A) = ipoo(A) = 2s(A) with 1 = 1,2. 
Further, using Euler-Maclaurin summation formula |l5], (13.261) becomes 

- / ^Ap(^)(A)pS;^(A) + — p(^)(A,)+ P^^'(A2)} (3.27) 

(I3.24P can also be expressed in similar form 

.Wfw , J2).,^ ^ ... r ^^Jl)f.^ 1 A A p'(A-Ai) 



P^^(A)+P)^^(A) - Poo(A)=/ rf7p(A-7)p);^(7)-TT77P(A-Ai) 



lAr ' ^ ^" 2iV 12iVVS^(Ai) 

+ />.p(A-o0.g'W-^p(A-A.)-^g^ (3,28) 

As before, pAi and pA2 are the largest sea roots from the two "seas" respectively. From 
this point, the calculation very closely resembles the details found in Section 2 in [3]. Hence, 
we omit the details and give only the crucial steps. Note that (13.281) can be written in the 
standard form of the Wiener-Hopf equation [16] after redefining the terms, 

/•OO /»oo 

X^'\t) + x^'\t) - / ds Pit - s)x^'\s) - ds Pit - s)x^'\s) 
Jo Jo 
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where the following definitions have been adopted 

x(')(A)=pS(a + aO 

/^'nA)=pS(A + AO (3.30) 

and following change in variable is used : t = X — Ai with I = 1 ,2 From the Fourier 
transformed version of fl3.29p . one can solve for X^\uj) which is the Fourier transfrom of 
{t) that is analytic in the upper half complex plane { 



2iV ■ 12iVV!;^(Az) 

igi 1 iuj 



+ 



Ll2iV2pff(A0 2iV 12iV2pJ)(A, 
Ufa 1 igi 



7r-tuj\N 2N i2Ar2pS(A, 



(3.31) 



where G+{uj)G+{-uj) = 1 + 02(0;), ^1 = ^^(2 + u - ^) and a = = (2(1!^) 2, with 

G+(0)^ = 

From (I3.3p . (13.301) and (I3.3ip . one can then determine p^''(Ai) and pj^''(A2) exphcitly from 

X? (0) = -p« {Ai) = - du; ) (.;) (3.32) 
by contour integration and some algebra. We give the result below 

pSi^(AO = ^jyT + 2™ + igi + [tt^ + - ^ + 4772^2 + 47ra(7r + z^^i)]^} (3.33) 

Finally, using pS(A) ~ 2e^'"^ for A — A; and (I3.27p . one arrives at the desired expression 
for 1/A^ correction to the energy, 

TT^ sin p 2\ 
En — Eoo = Ecasimir = — 24/I/V~ "~ ^'^'^ (3.34) 



where the effective central charge is 

"eff 



Cpff = 1 — 12a^ 



= 1-6^ (3.35) 

We see that for this model, the central charge, c = 1 (Free boson). Also Cg// is independent 
of boundary parameters, unlike for the Dirichlet case [3]. This is a feature expected for 



'Again for complete details, refer to [3] 
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models with Neumann boundary condition. Further, from conformal field theory, one also 
has the following for the conformal dimensions. 



A 



c, 



eff 



24 



4(z/-l) 



(3.36) 



Note that the above results are derived for the lowest energy state for even with one hole 
for each Qa{u). Reviewing the derivation above, one can notice that the results above can 
be further generalized for any and for low lying excited states with arbitrary number of 
holes, provided these holes are located to the right of the largest sea root as mentioned in 
the beginning of Section [31 For these excited states, the sum for S^^{X) in (13.231) - (I3.26P will 
inevitably have different limits since the number of sea roots vary. However, after applying 
the Euler-Maclaurin formula, one would recover (I3.27P and (13.281) . In addition to that, for 
states with number of holes (all located to the right of the largest sea root), one uses 
the more general result for the sum rule, namely (13.51) and (13.121) which eventually yields 

Nh 

Thus, we have the following for the effective central charge and conformal dimensions for 
low lying excited states 

c^rr = 1-6 , ^ 
^ii (z/ - 1) " 

Njj (3.38) 



4(^.-1)- 

Surprisingly, the results (13.360 and (13.38^ appear to have more resemblance to spin chains 
with diagonal boundary terms, as one could see from the dependance [33]- [36], rather 
than [12] which is the anticipated form for conformal dimensions for spin chains with 
nondiagonal boundary terms. Indeed the theory of a free Bosonic field (p compactified on a 
circle of radius r is invariant under ^ ip + 27rr, where r = |. (3 is the continuum bulk 
coupling constant that is related to u hy jS"^ = 87r(^^). Further, the quantization of the 
momentum zero- mode IIq, yields IIq = ^ for Neumann boundary condition and Hq = ^ for 
the Dirichlet case, where n is an integer. Hence, the zero-mode contribution to the energy, 
E^^n ~ Hq implies i?o,n ~ A ~ (^^) for Neumann and i?o,n ~ A ~ (^jz^) Dirichlet case 
respectively. More complete discussion on this topic can be found in [35l l42]. Next, we will 
resort to numerical analysis to confirm our analytical results obtained in this section. 
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4 Numerical results 



We present here some numerical results for both odd and even cases, to support our 
analytical derivations in Section [3^ We first solve numerically the Bethe equations (12. Sp . 
fl2.12p . fl2.13p . fl2.19p and (12.201) for some large number of spins. We use these solutions to 
calculate Casimir energy numerically from the following 

-^bulk ~l~ -^boundary ~l~ -^^Casimir (^•-^) 

In (14. ip . E is given by (13. ip . Thus, having determined the Bethe roots numerically, one uses 
known expressions for Ebuik [H] and Eboundary [38] to determine Ecasimir- Then using the 
expression found above for Ecasimir, namely (13.340 . one can determine the effective central 
charge, Cg// for that value of N, 

Ceff = TT^ {E — Ebulk — Eboundary) (4.2) 

Finally, we employ an algorithm due to Vanden Broeck and Schwartz [39]-[l0] to extrapolate 
these values for central charge at — > oo limit. Table 1 below shows the Cg// values for some 
finite even A^, for the lowest energy state with one hole (A'^^^ = 1). Equation (I3.38P predicts 
Ceff values of -11 and -7 for p = 1 and p = 3 respectively which are the extrapolated values 
(-11.000315 and -7.000410) we obtain from the Vanden Broeck and Schwartz method. 



A^ 


Ceff, p=l,iy = 2 


Ceff, p = 3,z/ = 4 


16 


-9.365620 


-2.853872 


24 


-9.857713 


-3.271279 


32 


-10.122128 


-3.557148 


40 


-10.287160 


-3.770882 


48 


-10.399970 


-3.939554 


56 


-10.481956 


-4.077652 


64 


-10.544233 


-4.193784 


oo 


-11.000315 


-7.000410 



Table 1: Central charge values, Cg// for p = 1 (a+ = 0.783, a_ = 0.859) and p = 3 
(a+ = 2.29, a_ = 1.76), from numerical computations based on A^ = 16,24,... ,64 and 
extrapolated values at A^ — > oo limit (Vanden Broeck and Schwartz algorithm). 

For odd A^ sector, since Nh = 0, (13.380 predicts Cg// = 1 (for the ground state) for any 
odd p. We present similar numerical results for odd A^ in Table 2 below for p = 1 and 

10^, 1 
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p — 3. We work out the Cg// values numerically for N — 15 ,25 ,. . . ,65. Excellent agreement 
between the calculated and the extrapolated values of 1.000770 and 1.001851 again strongly 
supports our analytical results. 



N 


Ce//, p=l,iy = 2 


Ceff, p = 3,U = i 


1 5 


n 8Q8334 




25 


0.936128 


0.634012 


35 


0.953433 


0.692758 


45 


0.963360 


0.731841 


55 


0.969797 


0.760142 


65 


0.974311 


0.781795 


oo 


1.000770 


1.001851 



Table 2: Central charge values, Cg// for p = 1 (a+ = 0.926, a_ = 0.654) and p = 3 
(a+ = 2.10, a_ = 1.80), from numerical computations based on = 15,25,... ,65 and 
extrapolated values at A'" ^ oo limit (Vanden Broeck and Schwartz algorithm). 

5 Discussion 

From the proposed Bethe ansatz equations for an open XXZ spin chain with special non- 
diagonal boundary terms at roots of unity, we computed finite size effect, namely the 1/N 
correction (Casimir energy) to the lowest energy state for both even and odd N. We also 
studied the bulk excitations due to holes. We found some peculiar results for these excita- 
tions of this model. Firstly, the number of holes for excited states seem to be reversed: even 
number of holes for chains with odd number of spins and vice versa. However, one could 
explain this by resorting to heuristic arguments involving effects of magnetic fields on the 
spins at the boundary. We then computed the energy due to hole-excitations. We further 
generalized the finite-size correction calculation to include multi-hole excited states, where 
these holes are situated to the right of the largest sea root. Having found the correction, we 
proceeded to compute the effective central charge, Cg// and the conformal dimensions, A for 
the model. We found the central charge, c — 1. The effective central charge is independent 
of the boundary parameters, as expected for models with Neumann boundary condition. 
The result for A however, turns out to be similar to models with diagonal boundary terms 
rather than the nondiagonal ones, to which the model studied here belongs to. 

As an independent check to our analytical results, we also solved the model numeri- 
cally for some large values of A^. We used this solution to compute 1/A^ correction for 
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these large N values, then extrapolate them to the — > oo limit using Vanden Broeck 
and Schwartz algorithm. Our numerical results strongly support the analytical derivations 
presented here. Spectral equivalences between diagonal-nondiagonal and diagonal-diagonal, 
nondiagonal-nondiagonal and diagonal-diagonal [T5[ |16[ |50] open XXZ spin chains have been 
shown to exist. Hence, one may attempt to explain the diagonal (Dirichlet) behaviour of the 
model studied here by some such equivalence. However, to our knowledge, such equivalences 
have been found when the boundary parameters obey certain constraint which is not 

the case for the model we considered here, as already remarked in Footnote 1. Hence, the 
question about the "Dirichlet-like" behaviour remains for now. We hope to be able to resolve 
this issue soon. 

There are many other open questions that one can explore and address further. For 
example, similar analysis involving boundary excitations can also be carried out. This can 
be really challenging even for the diagonal (Dirichlet) case [Mf HQ]. Further, solution for 
more general XXZ model involving multiple Q{u) functions [121 [13], can also be utilized in 
similar capacity to explore these effects. Last but not least, excitations due to other objects 
that we choose to ignore here, such as special roots/holes and so forth can also be explored 
for these models in order to make the study more complete. We look forward to address 
some of these issues in near future. 
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